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Abstract 

We report a search for narrow resonances, produced in pp collisions at ^fs = 1.96 TeV, that decay into muon pairs with 
invariant mass between 6.3 and 9.0 GeV/c 2 . The data, collected with the CDF II detector at the Fermilab Tevatron collider, 
correspond to an integrated luminosity of 630 pb _1 . We use the dimuon invariant mass distribution to set 90% upper credible 
limits of about 1% to the ratio of the production cross section times muonic branching fraction of possible narrow resonances 
to that of the T(1S) meson. 

PACS numbers: 13.20.Gd, 13.60.Lc 
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I. INTRODUCTION 



Searches for narrow resonances at SPEAR have set upper bounds on Ti, the leptonic width of possible resonances, 
of 100 eV in the cms energy region 5.7 < E cm < 6.4 GeV and of approximately 60 eV in the region 7.0 < E cm < 7.4 
GeV These limits do not exclude bound states of charge —1/3 and spin light quarks, the existence of which 
is required by supersymmetric theories. The leptonic width of £ mesons, 1 bound states of charge —1/3 scalar 
quarks, has been evaluated in Ref. Q using potential models of ordinary heavy quarks Q. Because of the p-wave 
suppression of the fermion contribution, the leptonic width is estimated to be approximately 18 (6) eV for a resonance 
with a 6 (10) GeV/c 2 mass. A more recent study [|| has set an average 90% upper credible limit of T/ ~ 8 eV to the 
existence of narrow resonances in the invariant mass range 6.3 — 9GeV/c 2 by using a dimuon sample corresponding 
to a luminosity of approximately 110 pb _1 collected with the CDF detector during the 1992— 1995 Tevatron collider 
run. Because of statistical fluctuations, that study does not rule out the existence of £ states over the full mass range. 
In particular, at a mass of 7.2GeV/c 2 , that study returns an upper credible limit that is consistent with what is 
expected if an £ state is produced. 

In this study, we repeat the analysis in Ref. [4] and take advantage of a higher statistics data sample that corresponds 
to an integrated luminosity of 630 pb _1 collected with the CDF II detector after May 2006. At that time, the CDF II 
trigger system was upgraded and made capable of acquiring events containing muon pairs with invariant mass larger 
than 6 GeV/c 2 and a kinematic acceptance comparable to that of the Run I trigger system. Section |TT] describes the 
detector systems relevant to this analysis. The data sample and analysis method are described in Sec. IIIII In Sec. IIV| 
we fit the dimuon invariant mass distribution and derive 90% upper credible limits to Ti as a function of the resonance 
mass. Our conclusions are summarized in Sec. El 

II. CDF II DETECTOR AND TRIGGER SYSTEM 

CDF II is a multipurpose detector, equipped with a charged particle spectrometer and a finely segmented calorime- 
ter. In this section, we briefly describe the detector components that are relevant to this analysis. A detailed 
description of these subsystems can be found in Refs. [j| @, 0, H, 0, [HI EH, E2, E3, EiJ- Two devices inside the 1.4 T 
solenoid are used for measuring the momentum of charged particles: the silicon vertex detector (SVXII and ISL) and 
the central tracking chamber (COT). The SVXII detector consists of microstrip sensors arranged in six cylindrical 
shells with radii between 1.5 and 10.6 cm, and with a total z coverage [l5| of 90 cm. The first SVXII layer, also 
referred to as the LOO detector, is made of single-sided sensors mounted on the beryllium beam pipe. The remaining 
five SVXII layers are made of double-sided sensors and are divided into three contiguous five-layer sections along the 
beam direction z. The two additional silicon layers of the ISL help to link tracks in the COT to hits in the SVXII. 
The COT is a cylindrical drift chamber containing 96 sense wire layers grouped into eight alternating superlayers of 
axial and stereo wires. Its active volume covers \z\ < 155 cm and 40 to 140 cm in radius. The transverse momentum 
resolution of tracks reconstructed using COT hits is ct(pt)Ipt — 0.0017 [GeV/c] -1 . 

The central muon detector (CMU) is located around the central electromagnetic and hadronic calorimeters, which 
have a thickness of 5.5 interaction lengths at normal incidence. The CMU detector covers a nominal pseudorapidity 
range \r)\ < 0.63 relative to the center of the detector, and is segmented into two barrels of 24 modules, each covering 
15° in (f>. Every module is further segmented into three submodules, each covering 4.2° in </> and consisting of four 
layers of drift chambers. The smallest drift unit, called a stack, covers a 1.2° angle in <fi. Adjacent pairs of stacks 
are combined together into a tower. A track segment (hits in two out of four layers of a stack) detected in a tower is 
referred to as a CMU stub. A second set of muon drift chambers (CMP) is located behind an additional steel absorber 
of 3.3 interaction lengths. The chambers are 640 cm long and are arranged axially to form a box around the central 
detector. The CMP detector covers a nominal pseudorapidity range \t]\ < 0.54 relative to the center of the detector. 
Muons that produce a stub in both CMU and CMP systems are called CMUP muons. The CMX muon detector 
consists of eight drift chamber layers and scintillation counters positioned behind the hadron calorimeter. The CMX 
detector extends the muon coverage to \r]\ < 1. 

The luminosity is measured using gaseous Cherenkov counters (CLC) that monitor the rate of inelastic pp collisions. 
The inelastic pp cross section at ^/s = 1960 GeV is scaled from measurements at y/s = 1800 GeV using the calculations 
in Ref. |l6[. The integrated luminosity is determined with a 6% systematic uncertainty [17| . 

CDF uses a three- level trigger system. At Level 1 (LI), data from every beam crossing are stored in a pipeline 
capable of buffering data from 42 beam crossings. The LI trigger either rejects events or copies them into one of the 
four Level 2 (L2) buffers. Events that pass the LI and L2 selection criteria are sent to the Level 3 (L3) trigger, a 
cluster of computers running speed-optimized reconstruction code. For this study, we select events with two muon 
candidates identified by the LI and L2 triggers. The LI trigger uses tracks with px > 1-5 GeV/c found by a fast track 
processor (XFT). The XFT examines COT hits from the four axial superlayers and uses the r — (j> information to find 
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tracks with pt > 1.5 GeV/c in azimuthal sections of 1.25°. The XFT passes the tracks to a set of extrapolation units 
that determine the CMU (CMX) towers in which a CMU (CMX) stub should be found if the track is a muon. If a 
stub is found, a LI CMU (CMX) primitive is generated. The LI dimuon trigger requires one CMU primitive and an 
additional CMU or CMX primitive. The L2 trigger additionally requires that at least one of the muons has a CMUP 
stub matched to an XFT track with pr > 3 GeV/c and that additional CMU (CMX) stubs are matched to XFT 
tracks with pt > 1.5 (2) GeV/c. The L2 trigger is dynamically prescaled as a function of the Tevatron luminosity 
in order to optimize the data acquisition performance. The L3 trigger selects muon pairs with invariant mass larger 
than 5 GeV/c 2 , and \5zq\ < 5 cm, where zq is the z coordinate of the muon track at its point of closest approach to 
the beam line in the r — <p plane. These requirements define the dimuon trigger used in this analysis. 

III. DIMUON DATA SAMPLE 

In this analysis, we select events acquired with the dimuon trigger and which contain two and only two muons 
with opposite charge. Events are reconstructed offline taking advantage of more refined calibration constants and 
reconstruction algorithms. COT tracks are extrapolated into the SVXII detector, and refitted adding silicon hits 
consistent with the track extrapolation into the microvertex detector. We require that at least one muon has stubs 
reconstructed in both the CMU and CMP detectors and is matched to a track with p? > 3 GeV/c. A track is 
identified as a CMUP muon if Arcj), the distance in the r — <f> plane between the track projected to the CMU (CMP) 
chambers and a CMU (CMP) stub, is less than 30 (40) cm. The additional muon is accepted if a stub reconstructed 
in the CMU or CMX detectors is matched to a track with p T > 2 GeV/c. A track is identified as a CMU (CMX) 
muon if Arcj) is less than 30 (50) cm. We require that muon-candidate stubs correspond to LI primitives, and correct 
the muon momentum for energy losses in the detector. 

Previous measurements of the B hadron mixing probability [l8| and of the bb correlated production [l9| have used 
similarly selected dimuon data samples. Those studies show that 70 — 75% of these dimuon pairs are contributed 
by bb and cc production. Following the procedure pioneered in Ref. Q, we reject muons arising from the decay of 
hadrons with heavy flavor with the following requirements: 

1. Isolation. The isolation, /, is defined as the scalar sum of the transverse momenta of all the tracks in a cone 
of radius R = y Scj) 2 + 6rj 2 = 0.4 around the muon direction. We require that both muons have isolation 
I < 4 GeV/c. 

2. Promptness. In contrast to b and c hadrons, the T(1S) and £ mesons have a negligible lifetime. We evaluate the 
invariant mass of each muon pair by constraining the two muon tracks to originate from a common point in the 
three-dimensional space (vertex constraint). We reject muon pairs if the probability of the vertex-constrained 
fit is smaller than 0.001. We also reject muon pairs if L xy /<j[ JT > 3, where L xy is the displacement of the 
dimuon-candidate vertex with respect to the pp event vertex projected onto the dimuon transverse momentum 
vector, and ax is its uncertainty. 

The invariant mass distribution of the muon pairs passing our selection requirements is shown in Fig. [T] 

IV. SEARCH FOR NARROW RESONANCES 

As in Ref. Q, we limit our search for narrow resonances to the mass region above 6.3 GeV/c 2 , where the kinematic 
acceptance becomes independent of the ^-candidate transverse momentum. The number of T mesons in the data 
is derived by fitting the invariant mass distribution in Fig. [T]with a binned maximum likelihood method [20]. We 
use a fifth-order polynomial to model the continuum in the invariant mass region 6 — 12 GeV/c 2 and Gaussian 
functions to model the T contributions l . The best fit returns 52700 ± 350 T(1S) meson events. The fit also returns 
-^T(iS) = 9459 ± 1 MeV/c 2 and a mass resolution au — 52 ± 1 MeV/c 2 . This mass resolution is well modeled by a 
simulation of the process pp — > T(1S)A". In the simulation, the event generator produces unpolarized T(1S) mesons 
with the transverse momentum distribution of the data [2l| and a flat rapidity distribution for |y| < 1. The generated 
events are processed with the CDF II detector simulation that in turn is based on the GEANT simulation package [13] • 
Events are then required to pass the same trigger requirements, analysis selection, and reconstruction criteria imposed 



Two Gaussian functions are required to model each T(1S) or T(2S) signal, whereas a single Gaussian function is sufficient to describe 
the T(3S) signal. 



6 



20000 




on the data. The same simulation 2 predicts an invariant mass resolution that increases from 32 MeV/c 2 for an £ 
state with a mass of 6.3 GeV/c 2 to 50 MeV/c 2 for m £ = 9 GeV/c 2 . 

The invariant mass distribution in Fig. [Tj does not exhibit any significant narrow bump in the mass range 6.3 — 
9 GeV/c 2 . We set limits to R = [erg BR(£ — > M/ i )]/[ cr T(is) Si?(T(lS) — > /u/i)], the production cross section times the 
muonic branching fraction, under the assumptions that £ states have negligible widths with respect to the experimental 
resolution and are produced unpolarized. We also assume that the efficiency of the isolation requirement, used to 
suppress the heavy flavor contribution, does not depend on the £ mass and is equal to that for T(1S) mesons. Limits 
on the existence of narrow £ states are set by adding a Gaussian term to the likelihood function used to fit the data in 
the 6—12 GeV/c 2 mass interval. We perform 108 fits in which we change the position of the Gaussian peak in steps 
of 25 MeV/c 2 from 6.3 to 9.0 GeV/c 2 . In each fit, we force the Gaussian width to the rms mass resolution predicted 
by the detector simulation for that mass. The usage of only one Gaussian function to estimate the number of £ events 
introduces a 4.5% systematic uncertainty. The uncertainty is estimated as the difference between the number of T(1S) 
events returned by fits that use one or two Gaussian functions to estimate the T(1S) signal. For each £ mass value, 
we use the integral of the Gaussian function and its error returned by the best fit to derive N u i, the 90% credibility 
upper limit to the number of events contributed by a narrow resonance 3 . We evaluate the ratio of the geometric 
and kinematic acceptance for an £ resonance to that for the T(1S) meson with the previously described Monte Carlo 
simulation. The geometric and kinematic acceptance depends on the resonance mass and increases from 65.5% of that 
of the T(1S) meson for mg- = 6.3 GeV/c 2 to 97.4% for mg = 9.0 GeV/c 2 . Small data-to-simulation correction factors 
for the trigger and reconstruction efficiency p^ . |23| depend little on the muon pt and are neglected as in the study 
in Ref. [4J. The kinematic acceptance for T(1S) mesons depends on their polarization. If T(1S) mesons are produced 
longitudinally polarized, the acceptance increases by 5%, whereas it decreases by 3% if T(1S) mesons are produced 
transversely polarized. The current knowledge of the polarization of T(1S) mesons produced at the Tevatron [241. |25| 
yields an additional 3% systematic uncertainty on the ratio of the £ to T(1S) production cross sections. The total 
systematic uncertainty on the ratio (6%) is added in quadrature to the statistical errors returned by the fits. 



2 When simulating £ states, the transverse momentum distribution is rescaled from that of the T(1S) data so that < > / < 
~f( IS) 

p T >= mg-/m-f (is) ■ For J/ip mesons, this rescaling procedure predicts a da/dpx distribution that decreases more rapidly with 
increasing momenta than the distribution of the data. However, a poor modeling of the transverse momentum distribution is not a 
cause of significant error because the kinematic acceptance does not depend on the £ transverse momentum within a couple of percent. 

3 The integral of the fit likelihood from N u i to infinity is 10% of the integral of the fit likelihood from to infinity. The justification for 
this procedure is Bayesian with a prior that is zero for negative resonance cross sections and flat for positive ones. 
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FIG. 2: Bayesian 90% upper credible limits (histogram) to a ^ ,-, ,, Ti|< 

dashed histogram represents limits that include the 6% systematic uncertainty. The solid line represents the limits expected 
when fits return no £ events. 
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FIG. 3: Bayesian 90% upper limits to Tf (histogram). The dashed histogram represents limits that include the 6% systematic 
uncertainty. The solid line represents the limits expected when fits return no £ events. The dashed and dotted lines represent 
the leptonic widths of I bound states of scalar quarks predicted in Refs. Q] and |26|. respectively. 



The 90% credibility upper limit to R is given by the ratio of N u i to the number of observed T(1S) mesons, 
corrected for the relative acceptance. These limits are shown in Fig. [5] as a function of the £ mass. Figure [3] shows 
the 90% upper credible limits to Fp. The latter limits are derived under the assumption 4 that <J£ BR{£ — > /x/x) = 



4 As noted in Ref. Q|, this equation tends to underestimate the production cross section of states lighter than the T(1S) meson. If this 
were the case, the r? limits set by our study would be correspondingly smaller than those indicated in Fig. [3] 
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<7 T( is) BR(T(1S) -> MM) x (m T(ls) /m £ ) 3 x rf /rJ (1S) 0. 



V. CONCLUSIONS 



We present a study of the invariant mass spectrum of muon pairs produced at the Fermilab Tevatron collider and 
collected by the CDF experiment. This study improves previous limits to the existence of narrow resonances in the 
invariant mass range 6.3 — 9.0 GeV/c 2 . We derive 90% upper credible limits to R, the ratio of the production cross 
section times muonic branching fraction of possible narrow resonances to that of the T(1S) meson, of the order of 
1%. Assuming that erg BR(£ ) = or(is) 

RR(T(1S) -> mi) x (m T[ls) /m £ ) 3 x T^/T^ 10 ', these limits correspond to 
90% upper credible limits smaller than 10 eV to the leptonic width of possible narrow resonances. 
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